Available online at www.sciencedirect.com 


ScienceDirect 


JOURNAL OF 


www.elsevier.com/locate /jpowsour 


Journal of Power Sources 163 (2007) 907-915 


A three-dimensional two-phase flow model for a 
liquid-fed direct methanol fuel cell 


Jiabin Ge, Hongtan Liu * 


Department of Mechanical and Aerospace Engineering, University of Miami, Coral Gables, USA 


Received 10 August 2006; received in revised form 3 October 2006; accepted 4 October 2006 
Available online 13 November 2006 


Abstract 


A three-dimensional, two-phase, multi-component model has been developed for a liquid-fed DMFC. The modeling domain consists of the 
membrane, two catalyst layers, two diffusion layers, and two channels. Both liquid and gas phases are considered in the entire anode, including the 
channel, the diffusion layer and the catalyst layer; while at the cathode, two phases are considered in the gas diffusion layer and the catalyst layer 
but only single gas phase is considered in the channels. For electrochemical kinetics, the Tafel equation incorporating the effects of two phases is 
used at both the cathode and anode sides. At the anode side the presence of gas phase reduces the active catalyst areas, while at the cathode side 
the presence of liquid water reduces the active catalyst areas. The mixed potential effects due to methanol crossover are also included in the model. 
The results from the two-phase flow mode fit the experimental results better than those from the single-phase model. The modeling results show 
that the single-phase models over-predict methanol crossover. The modeling results also show that the porosity of the anode diffusion layer plays 
an important role in the DMFC performance. With low diffusion layer porosity, the produced carbon dioxide cannot be removed effectively from 
the catalyst layer, thus reducing the active catalyst area as well as blocking methanol from reaching the reaction zone. A similar effect exits in the 


cathode for the liquid water. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Efforts in developing mathematical models for DMFC have 
been limited until recent years. Baxter et al. [1] developed a one- 
dimensional, single-phase mathematical model for a liquid-fed 
DMFC, focused on the anode catalyst layer. Dohle et al. [2] 
presented a one-dimensional model for the vapor-fed DMFC 
including methanol crossover and the effects of methanol con- 
centration on the cell performance were studied. Scott et al. 
[3,4] developed several simplified single-phase models to study 
transport and electrochemical processes in the liquid-fed DMFC 
and showed that the cell performance was limited by the slow 
diffusion of liquid methanol. Sundmacher and Scott [5] devel- 
oped a model of the methanol mass transport process and the 
model was used to predict the effective methanol concentration 
at the catalyst surface and the anode polarization. This model, 
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together with an empirical model of the open circuit voltage 
and the cathode overpotential model, was used to predict the 
overall cell voltage and current density response of the fuel 
cell. Cruickshank and Scott [6] presented a simple model to 
describe the permeation of methanol from anode to cathode in 
DMEC. Shukla et al. [7] developed a one-dimensional model for 
a comparison of the performance of a solid—polymer electrolyte 
DMFC with aqueous methanol and methanol mixed with air at 
the anode. Kulikovsky [8] simulated a vapor-fed DMFC with a 
two-dimensional model and compared the detailed current den- 
sity distributions in the backing, catalyst layer, and membrane 
separator between a conventional and a new current collector. 
Wang and Wang [9] developed a two-dimensional, two-phase 
and multi-component model for a liquid-feed DMFC. In the 
model, the catalyst layer was considered as an infinitely thin 
interface. Birgersson et al. [10] developed a two-dimensional 
single-phase DMFC anode model, where the catalyst layer was 
considered as boundary conditions via parameter adaptation. 
Divisek et al. [11] accounted for the influence of both the 
methanol and oxygen reaction kinetics and their dependency 
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Nomenclature 


C mass fraction 

D diffusivity (m? s7!) 

F Faraday constant (96485.309 C mol!) 

k permeability of gas diffuser (m7) 

M molar mass (kg mol`!) 

P pressure (Pa) 

p% the saturation pressure of water at operating tem- 


perature 
r the porous media correction factor 
Te advection correction factor 
R gas constant (8.314 J (mol K)"!) 
Rmr resistance of membrane (Q m) 
sS saturation 
S source term 
V the intrinsic velocity vector (m s7!) 
X mole fraction 
z charge of the fixed sites 


Greek letters 


porosity 
0) membrane phase potential 
À the electro-osmotic drag coefficient 
H viscosity (kg (m s)7!) 
v kinetics viscosity (m? s7!) 
Oe contact angle 
p density (kg m~?) 
o surface tension (N m~ ty 
E relative mobility 


Subscripts and superscripts 
anode, air 
cathode 
diffusion 
effective 

gas 

proton 

ith layer 

kth component 
liquid water 
methanol 
membrane 
nitrogen 
oxygen 

water 


Td CAO & 


+ 


£O2B BTS 


on the two-phase mass flow in the catalyst and diffusion layers 
by a vapor-—liquid model for the DMFC. Detailed information 
on fuel cell modeling and transport phenomena can be found in 
several recent reviews [12—14]. 

It is well-known that reaction rate in a catalyst layer is not 
uniform. Therefore, to accurately modeling a fuel cell, including 
a DMFC, catalyst layers must be treated as parts of the model- 
ing domain instead of boundaries or interfaces. Furthermore, a 
DMFC is a complex multiphase and multi-component system. 


WC 
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Fig. 1. The schematic of the modeling domain. 


Multiphase flows are driven by buoyancy, gravitational, capil- 
lary and viscous forces. For multiphase flows in porous media, 
the capillary forces play fundamental roles in controlling phase 
distributions. A model for a DMFC without including the two- 
phase flow is not complete. 

In this work, a three-dimensional two-phase model is devel- 
oped and the modeling domain includes the two fluid channels, 
two diffusion layers and two catalyst layers. The model is used 
to study various transport phenomena in a DMFC, the cell per- 
formance and the effects of methanol crossover. 


2. Mathematical model development 


Fig. 1 shows the modeling domain for a DMFC, which con- 
sists of the anode side, cathode side and polymer electrolyte 
membrane. The anode side and cathode side have the same struc- 
ture, which is divided into the gas channel, diffusion layer and 
catalyst layer. The fuel and oxidant flow along the channels, 
where the flows are assumed to be laminar and the gas mixtures 
are assumed to be perfect gases. All phases are assumed to be 
continuous so that the continuum approach can be applied. In 
the cathode side, two phases are considered in the gas diffusion 
layer and the catalyst layer, but in the gas channels, only gas 
mixture of water vapor, oxygen and nitrogen is considered and 
the liquid water is neglected. In the anode side, two phases, the 
liquid phase of mixture of water and methanol, and gas phase 
of carbon dioxide are considered in all the areas, including the 
channels, the diffusion layer and the catalyst layer. 

Such parameters as gas volume fraction and liquid saturation 
are assumed to be volume-averaged properties to accommodate 
a macroscopic continuum approach in the porous media. Liq- 
uid water is assumed to be homogeneously distributed within 
each control volume. The liquid phase velocity is different from 
the gas phase velocity and the mixture velocity is the weighted 
average of these two. 

In Fig. 1, hc is the gas channel height; td, tc, tm are the thick- 
ness of the diffusion layer, catalyst layer, and the membrane, 
respectively; wc is the half-width of the gas channel; ws is the 
width of the collector plate; and Lc is the length of the fuel 
cell. The reactants of anode side and cathode side enter the gas 
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Table 1 

Geometric parameters used in the model 

Channel length, Lc (m) 6.5 x 107? 
Half-channel width, we (m) 5.0 x 1074 
Channel height, hc (m) 8 x 1074 
Channel shoulder width, ws (m) 1.0 x 1073 
Diffusion layer thickness, td (m) 1.4 x 1074 
Catalyst layer thickness, tc (m) 3 x 1075 
Membrane thickness (Nafion® 117), tm (m) 1.8 x 1074 


channels from the surface at y=0. The geometric parameters are 
given in Table 1. 


3. Governing equations 


The governing equations include mass conservation equa- 
tion, momentum equations, species conservation equations, and 
electrochemical kinetics equations (Tafel equation). All the con- 
servation equations are developed for the entire solution domain 
of the model with source terms modified in each region to reflect 
the appropriate physical phenomena. These equations are simi- 
lar to those used in PEM fuel cells [15,16], 


alep) 


V -(epV) =0 1 
ag TV EV) (1) 
a(epV 55 3 
7 ) +V-(e0VV) = -EVP+V-(euVepV)+ pS; D 
V -(yepVC) = V - (psi Di VC) + V - [Pgs DV Cg] 


=V- (C1 — Cg) Si] + Sx (3) 


where p is the total mixture density of liquid and gas phase, v 
the intrinsic velocity vector, ¢ the porosities, P the pressure, Sx 
the source terms, jz viscosity, sı liquid saturation and sg is the 
saturation of gaseous phase (void fraction). Note that saturation 
is defined as the volume fraction of a phase in the mixture and 
S1 +Sg= 1. 

Eq. (1) is the mass conservation equation for the entire cell. 
Eq. (2) is the momentum equation in vector form, where the 
source term accounts for the forces exerted on the fluid by the 
solid matrix of a porous media and is zero in the channel. Eq. 
(3) is the species conservation equation, where the source term 
accounts for the species created/consumed by electrochemical 
reaction in catalyst layers and is non-zero in catalyst layers only. 
In the cathode diffusion layer and catalyst layer, species equa- 
tions are given for oxygen, nitrogen, water vapor and liquid 
water; while in the anode side, species equations are given for 
liquid water, liquid methanol, dissolved carbon dioxide and car- 
bon dioxide gas. 

The density of the mixture of liquid and gas phase is given 
by, 


P = PIS] + PgSg (4) 
The water saturation is given by, 


E Pe(C™ — Coq) 
PICe — C") + Oel” = CY oq) 


(5) 


Sw 


and when 


CY < CY 


sa Sw =0 


where CY 


z eq 18 the equilibrium water concentration and it is given 
by, 


WwW 

cy, = ROD m 

' Pg RT 
where P,(T) is the water vapor saturation pressure at the given 
local temperature. 

In Eq. (3) Dy (k=1or g) is the effective diffusivity of compo- 
nent kin the porous media and it is related to its binary diffusivity 
by the following equation: 


DÈ = Dye! (7) 
The advection correction factor is given by, 


— PECI + EgCg) 


= (8) 
° pisıCı + PgSgCg 
and the relative mobility are given by, 
ka/vı 
&(s)) = (9) 
(ka/™i) =P (krg/Vg) 
Eg(s1) = 1 — &(s1) (10) 
The relative permeabilities for liquid and gas phase are, 
ky = 8 (11) 
kre ma? 12 
Tg — Sg ( ) 
The fluid viscosity is given by, 
p 
u= (13) 
(ka/vi) + (krg/Ug) 
J; in Eq. (3) is given by, 
ElÉgk > 
J = E [V Pe + (01 — Pg)8] (14) 


where the capillary pressure saturation function, Pe, depends on 
the interfacial tension between the liquid water and the solid 
phase of the porous media. It can be determined by the pore 
geometry [15,17], 


e\1/2 
Po(s) = Py — A = o coste (=) Ys) (15) 


where ø is the surface tension at the gas liquid interface, the con- 
tact angle 6, of a liquid droplet defined as cos 6. =0sg — og//o, 
Pg the pressure of gas phase, P| the pressure of the liquid phase, 
and y(s) is the capillary pressure function. 

The cubic Leverett function is usually used for the capillary 
pressure function, 


W(s) = 1417s, — 2.12055 + 1.2635; (16) 


910 J. Ge, H. Liu / Journal of Power Sources 163 (2007) 907-915 


Table 2 

The source terms in the momentum equation 

Channel Te, 1.0 Si 0 

1 d bo 

Diffuser 2.25 — — 1 -eag V 
1 2 Hz ke 

Catalyst 2.25 ——1 —£c— V + —Cp+zFVo 
Ec k kh 


3.1. Reduction of the governing equations 


For steady state, the mass conservation and momentum equa- 
tions are reduced to, 


V -(epV) =0 (17) 


pV -VV =-VP4+rV (UVV) + pS; (18) 


where the source terms are given in Table 2. 

In the momentum equation and Table 2, ¢q and & are the 
porosities of the diffusion layer and catalyst layer, respectively, 
rs is the porous media correction factor, k is the permeability of 
water (ky) in the anode and the permeability of air (kair) in the 
cathode, kg is electro-kinetic permeability in the membrane, kp 
is the hydraulic permeability in the membrane, z is the charge 
number of the fixed sites, Cy+ is the concentration of protons 
that is taken to be the concentration of the fixed charge sites, and 
¢ is the membrane phase potential. 

The individual phase velocities are calculated from, 


ep\Vi = J + &epV (19) 
EpgVg = — J + &ge0V (20) 


3.2. The species equation in the cathode side 


Two gas species, oxygen and nitrogen, and the mixture of 
vapor and liquid water, are considered in the cathode. 

Assuming only gas phase (oxygen, nitrogen and vapor) exists 
in the cathode gas channel, the species equation is in the same 
form to that for single-phase flow (Table 3), 


In the cathode gas diffusion layer and catalyst layer, two-phase 
flow is considered and the species are oxygen (CO), nitrogen 
(CN) and water (C). 

When the dissolved oxygen in water is neglected, its concen- 
tration in the liquid phase is zero, then, 


Pe 


ce =0, 
PgSg 


pc? = psCo se and Fe = 


From Eq. (3), the oxygen species equation becomes, 


V - (Pg VgC9) = V - [PgsgD VC] + So (24) 


In the porous cathode electrode, the transport of liquid water 
is driven by two mechanisms, capillary force (driven by satura- 
tion gradient) and interfacial shear force (exerted by gas flow) 
[18-20]. The net water flux through the interface between the 
cathode catalyst layer and the GDL include the water generation 
rate in the catalyst layer and the net water transfer rate due to 
electro-osmosis and diffusion, 
w Mw(l +a). 


Pig 


where a is the net water transport coefficient across the mem- 
brane. Water exists in both vapor and liquid phases. For the liquid 
phase, the only species is water, so, 


cCY=1,  VC¥=0, 


(25) 


and pC = psy + PgSgCy 


The species equation for water becomes, 


V - (Pe VeCy) = V - [Pgs Dy VCZ] - V- (AV) + Sw 26 
c°+0N +C =1 (27) 


The source terms in the species equations are given in Table 3. 
3.3. The species equations in the anode side 


Carbon dioxide gas management is an important issue in 
liquid-fed direct methanol fuel cell [21]. The relatively large 
amount of carbon dioxide reduces the free area for methanol to 
reach the catalyst layer. Carbon dioxide is mostly transported out 
to the anode channel through the catalyst layer and the diffusion 


vy rO O 
V. (pg VgCg ) =V. (PgD9 VCS) (21) layer. 
-z For the liquid phase, the species include water, methanol and 
V - (Pg VgCg ) = V «(og Dg VCg ) (22) the dissolved carbon dioxide; while for the gas phase, only car- 
o m N bon dioxide is considered. 
Cg tCg tC =! (23) Thus, CS = 
Table 3 
Anode and cathode source term (s+) in the species equations 
Anode (source term sg) Cathode (source term sg) 
Methanol Water Proton Carbon dioxide Oxygen Water 
Channel layers 0 0 0 0 0 0 
Diffusion layers 0 0 0 0 0 0 
Catalyst | Mm . My . 1 Me , Mo My(1+a) . 
atalyst layers -j =A = eae aos. SWS 
y y 6F ” 6F 7* F” F” 4F” oF e 
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The liquid saturation is given by, 


1-—C 
ie Pst ) (28) 
C CE a) + Pell — C9) 
sı=1if C° < C? 


l,sat 
For methanol, Cy = 0, and the species equation for liquid 


methanol is, 
V- (AVC = V - (psiD§, VCP) + Sm (29) 


The mass fraction of carbon dioxide in the gas phase is C = 
1. 
The species equation for carbon dioxide is, 


V - (MCD = V - (psi DEVCf) — V - (pg Vg) + Se (30) 
The total mass fraction of carbon dioxide is, 

pC? = ps\Cy + pg(1 — sı) (31) 
And in the anode, we have 


cCo+C%+C™=1 (32) 
3.4, Electrochemical equations in catalyst layers 


Tafel equation is used to describe reaction in both the cath- 
ode and anode catalyst layers. In the anode catalyst layer, the 
carbon dioxide gas block methanol from reaching the catalyst 
surface and thus reduce the catalyst active surface area. In the 
cathode catalyst layer, the presence of liquid water contributes 
to electrode flooding by blocking active sites on the catalysts 
and reducing the gas volume in the porous electrode [22]. This 
leads to smaller active area available for reaction and reduced 
area and volume for oxygen gas transport [9,18]. To account for 
these coverage phenomena, the term sı is included in the Tafel 
equation for the anode side and the term (1 — sw) is included in 
the Tafel equation for the cathode side. Thus, 


i „ef 51C% Aa F 
ja = aif cref exp (25%) (33) 
c c 
(1 — sw)PgC? (e F ) 
. : -O 5g c 
Je + Jp = J0,ref exp n (34) 
c P re Pg,retC9 o RT c 


where jp is the pseudo-current accounting for methanol 
crossover flux. 


3.5. Cell potential and current density 


The cell voltage is calculated by, 
Ecen = Eo — na — Ne — [Ror (35) 


where Eo is the difference between the half cell potentials of the 
anode and cathode, Rmr the membrane resistance, and Z is the 
average current density of the cell and is given by, 


tc 
fe I re (36) 
0 


where tc is the catalyst layer thickness. 


The ionic conductivity of Nafion® is given by [23], 


Omel(T) = o'! exp | 1268 Ja (37) 
ce = 303 T 

where of! is the reference ionic conductivity at 303 K and is 

given by, 

ott = 0.005136 — 0.00326 (38) 


where the water content in the membrane, 6, depends on the 
water activity, a, according to the following relationship [24]: 


B = 0.043 + 17.18a — 39.85a? + 36.003, (0<a <1) 
B= 14.0+14@-1), (<a<3) 
The water activity is calculated by, 
a XYP (39) 
P(T) 


where X™ is the water mole fraction. 
The molecular weight of the mixture is given by, 


1 
M = =~ 40 
Em = 
The water mole fraction is given by, 
CYM 
XY = (41) 
MY 
The resistance of membrane is defined by, 
tm 1 
Rm = f — dz (42) 
0 Omr 


where tm is the membrane thickness. 
3.6. Methanol crossover 


In the membrane, the methanol transports due to diffusion 
and electro-osmotic drag and the flux is given by [7], 
mr __ Amt 1.5 mr dC 
m 7 F ~ m” m dz 
where Ch is the methanol concentration in the membrane. 
The electro-osmotic drag coefficient (Am) of methanol is 


defined as the number of methanol molecules dragged by each 
proton through membrane (CH30H/H*), and it is given by [25], 


Am = Xmlz=0Aw (44) 


(43) 


where àw is the electro-osmotic drag coefficient for water. 
The diffusion coefficient of methanol in membrane (D) is 
given by [7], 


D™ = 4.012 x 107" exp(0.0243127) [m* s~"] (45) 
3.7. Boundary conditions 


At the anode and cathode inlets, 


v= Vied? cr = Cied cy = Cred (anode) 
V? = Vied’ C7 = Co. Cy = Cieq (cathode) 
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Table 4 

The parameters and properties used in the model 

Property Value 

Cell temperature, T (K) 343 
Pressure, P (N m~?) 1.013 x 10° 
Anode reactant flow rate, ginlet (ml min™!) 15 

Cathode reactant flow rate, ginlet (ml min7!) 1000 

Inlet oxygen mole fraction, X ihlet (mol mol~!) 0.21 

Inlet methanol concentration, giaet (mol m~?) 1000 
Methanol reference concentration, es (mol m7?) 2000 
Oxygen reference concentration, ce (mol m~?) 0.472 
Diffusion coefficient of oxygen in air, Do (m? s7!) 1.22 x 107! 
Diffusion coefficient of water in air (m?) 3.32 x 10-5 
Diffusion coefficient of methanol in water, Dm (m? s7!) 2.8 x 107° 
Porosity of diffusion layer, eq 0.5, 
Porosity of catalyst layer, £c 0.6 
Porosity of membrane, £m 0.28 
Permeability to air in the gas diffuser, kair m?) 1.76 x 107!! 
Permeability to water in the gas diffuser, ky (m7) 1.0 x 107!! 
Air viscosity, Wg (kg (m s)7) 2.05 x 1075 
Electrokinetic permeability of the membrane, kg (m?) 7.18 x 107? 
Hydraulic permeability of the membrane, kp (m?) 1.8 x 10718 
Reference exchange current density times specific area at anode, airtel (Am?) 1.0 x 108 
Reference exchange current density times specific area at cathode, ait (A m~?) 200 

Anode reaction order, ya 1.0 
Cathode reaction order, ye 1.0 
Electro-osmotic drag coefficient of water, Ay, 25 

Anode transfer coefficient, œa 0.5 
Cathode transfer coefficient, œe 1.0 

Charge of fixed (sulfonate) sites, z —1 

Air density, Pair (kg m~?) 1.025 
Water density, Pwater (kg m7?) 977.8 

Net water transport coefficient, a 0.5 


At the walls and the interface with the collector plate shoulder, 
no-slip velocity boundary conditions are used. Neumann-type 
boundary conditions are used for velocity at the channel outlets. 
At the interface between the catalyst layer and membrane, the 
velocity is assumed to be zero. For all species equations, the 
Neumann boundary conditions are is applied at all the outlets, 
solid walls and symmetry surfaces. 


4. Modeling results and discussions 


4.1. Comparing with experimental data and single-phase 
model 


All parameters and properties in the model are given in 
Table 4. To validate the numerical model developed in the pre- 
ceding section, typical modeling results are compared with our 
experimental data [26] for a single cell operated at 1 M methanol 
concentration as shown in Fig. 2. The cell temperature is 70°C, 
cathode humidification temperature 70°C, methanol flow rate 
6 ml min“ !, and air flow rate 1200 SCCM. It is seen from Fig. 2 
that the result from the two-phase model shows better agreement 
with the experimental data than the single-phase model, partic- 
ularly at high current density region. Fig. 3 shows the modeling 
results of anode overpotential versus current density from the 
two phase model and the single-phase model at the same oper- 
ating parameters as those in Fig. 2. In the single-phase model, 


only the mixture of water and methanol is considered. Carbon 
dioxide is created by electrochemical reaction, partly dissolved 
in the liquid phase. When the dissolved carbon dioxide concen- 
tration exceeds its solubility, carbon dioxide gas is emerges and 
partially covers the reaction sites. From Fig. 3, it can be seen that 
at lower anode overpotential, the corresponding current density 


4 Experimental data 
Two-phase model 
Single-phase model 


0.8 


o 
D 


Voltage (V) 


0.1 0.2 0.3 0.4 0.5 0.6 
Current density (A/cm?) 


Fig. 2. Comparison of modeling results with experimental data and the results 
from single-phase model. Methanol feed concentration 1M, cell tempera- 
ture 70°C, cathode humidification temperature 70°C, methanol flow rate 
6ml min7~!, and air flow rate 1200 SCCM. 
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—t— Two-phase model 
—+:— Single-phase model 


Anode overpotential (V) 
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Current density (A/cm?) 


Fig. 3. Comparison of predicted of anode overpotential between the two-phase 


model and single-phase model. Methanol concentration 1 M, cell temperature 
70°C, methanol flow rate 6 ml min~!, and air flow rate 1200 SCCM. 
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Fig. 4. Predicted methanol crossover pseudo-current density as function of cell 
current density. Methanol concentration 1 M, cell temperature 70°C, methanol 
flow rate 6 ml min7!, and air flow rate 1200 SCCM. 


is low and the deviation of single-phase model result from that of 
the two phase model is negligible; while at high current density 
the single-phase model under-predicts the anode overpotential 
and the amount of under prediction increases with current den- 
sity. 


4.2. Methanol crossover 
Fig. 4 shows a modeling result for methanol crossover. 


The flux of methanol crossover is represented with a pseudo- 
current density. The methanol crossover is mainly caused by 


0.184 —ūO— Two-phase model 
i —*— Single-phase model 
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Fig. 5. The comparison of methanol crossover from the two-phase model 
and single-phase model. Methanol concentration 1 M, cell temperature, 70°C, 
methanol flow rate 6 ml min7!, and air flow rate 1200 SCCM. 


diffusion and electro-osmotic drag. Diffusion is caused by the 
methanol concentration gradient across the membrane. The 
electro-osmotic drag is caused by proton transfer through the 
membrane. When protons transfer from the anode to the cath- 
ode, water is dragged along. Since methanol is mixed with water 
and thus it is also dragged from the anode to the cathode. There- 
fore, drag coefficient for methanol is given by Eq. (44), from 
which it is clear that methanol flux due to electro-osmotic drag 
is proportional to its concentration at the anode/membrane inter- 
face. When the cell current density is about 0.073 Acm7?, the 
pseudo-crossover current density reaches its maximum. This 
phenomenon is caused by the dual effect of current density on 
methanol crossover. The diffusion flux decreases with current 
density due to the consumption of methanol in the anode catalyst 
layer, causing a lower methanol gradient across the membrane. 
The flux due to electro-osmotic drag is proportional to both the 
concentration at the anode and the cell current density. As current 
density increases, the product of methanol concentration and Aw 
may reach a maximum. Fig. 5 shows an example of comparison 
between the results of methanol crossover pseudo-current den- 
sities from the single-phase model and the two-phase model. It 
is clear that the single-phase model significantly under-predicts 
the methanol crossover. 


4.3. Anode liquid saturation 


Fig. 6 shows a contour plot of anode liquid saturation at cur- 
rent density of 0.367 A cm7? and anode overpotential of 0.52 V. 
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Fig. 6. Contour plot of anode liquid saturation. Cell current density, 0.367 A cm~? and anode overpotential 0.52 V, methanol concentration, 1 M, and cell temperature, 


70°C, methanol flow rate 6 ml min~!, and air flow rate 1200 SCCM. 
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Fig. 7. Contour plot of anode liquid saturation. Cell current density, 0.192 A cm7? and anode overpotential 0.4 V, methanol concentration, 1 M, and cell temperature, 


70°C, methanol flow rate 6 ml min~!, and air flow rate 1200 SCCM. 
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Fig. 8. The contour plot of total carbon dioxide mass fraction at current density 0.367 A cm~’, anode overpotential 0.52 V, and anode GDL porosity of 0.6, methanol 


flow rate 6 ml min~!, and air flow rate 1200 SCCM. 
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Fig. 9. The contour plot of total carbon dioxide mass fraction at current density 0.51 A cm, anode overpotential 0.52 V, the anode GDL porosity of 0.7, methanol 


flow rate 6 ml min~!, and air flow rate 1200 SCCM. 


Since the gas saturation satisfies sı + sg = 1, it is clear that gas 
volume fraction increases along the channel due to the accu- 
mulation from upstream. The modeling result also showed that 
the gas volume fraction in the catalyst layer can be as high as 
more than 40%. With such a high gas volume, a significant frac- 
tion of the catalyst surface can be made inactive. Fig. 7 shows 
another contour plot of anode liquid saturation at current density 
of 0.192 A cm7? and anode overpotential 0.40 V. Comparing it 
with Fig. 6, it can be seen that higher liquid fraction or lower 
carbon dioxide gas fraction exists at lower current density. 


4.4. Carbon dioxide mass fraction 


Fig. 8 shows the distribution of carbon dioxide mass fraction 
along the channel direction. Note that the mass fraction includes 
the masses of both the dissolved and gaseous carbon dioxide. The 
carbon dioxide mass fraction at the inlet is zero, and along the 
channel as electrochemical reaction occurs and the methanol 
is consumed, the carbon dioxide mass fraction increases. The 
amount of carbon dioxide created is proportional to the current 
density. Firstly, carbon dioxide dissolves in water and when it 
reaches its solubility limit, carbon dioxide gas emerges. At the 
same flow rate, the higher the current density is, the higher is 
the carbon dioxide mass fraction. 

Fig. 9 shows another result of distribution of carbon dioxide 
mass fraction. All the parameters are the same as those in Fig. 8 
except the porosity of the diffusion layer, which is 0.7, greater 
than that in Fig. 8. Due to the higher porosity of the diffusion 
layer, the corresponding current density is higher. Comparing 


with Fig. 8, it can be seen that the case shown in Fig. 9 has 
lower carbon dioxide concentration in the catalyst layer even 
though it has a higher current density. This means that a higher 
porosity in the diffusion layer lead to more effective removal of 
carbon dioxide gas, thus reducing the blockage of the catalyst 
surface area, resulting in higher current density and better cell 
performance. 

Fig. 10 shows the distribution of CO2 mass fraction along z- 
direction and at the different locations along the flow direction. 
It can be seen that near the channel inlet, the CO2 mass fraction 
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Fig. 10. Distribution of CO2 mass fraction along z-direction at different locations 
along the flow direction. Methanol concentration | M, cell temperature, 70°C, 
methanol flow rate 6 ml min7!, and air flow rate 1200 SCCM. 
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Fig. 11. Average CO2 mass fraction profiles at different current densities. 
Methanol concentration 1M, cell temperature, 70°C., methanol flow rate 
6 ml min7!, and air flow rate 1200 SCCM. 


is higher in the catalyst layer due to the higher current density. 
Since current density decrease along the channel direction the 
CO2 mass fraction also decreases along the channel direction. 
But CO2 mass fraction in the channel is lower near the inlet 
and increase along the flow direction. This is caused by the 
accumulation of the generated CO2. Fig. 11 shows the average 
CO2 mass fraction variation along the z-direction at different 
average cell current densities. The averaged values are taken 
over the entire channel, form the inlet to the outlet. It is clear 
that the higher the current density is the higher is the CO2 mass 
fraction. Since CO3 is generated in the catalyst layer and transfer 
out through the diffusion layer to the channel, the highest CO2 
mass fraction occurs in the catalyst layer. The dominant mass 
transfer mode in the diffusion layer is diffusion, so the profile 
of CO2 mass fraction in the diffusion layer is almost linear even 
after taking the average value. 


5. Concluding remarks 


A three-dimensional, two-phase, multi-component model has 
been developed for a liquid-fed DMFC. The anode and cathode 
catalyst layer are considered as parts of the modeling domain 
instead of interfaces. The model incorporates the effects of the 
second phase on the reduction of active catalyst surface areas 
and the mixed potential effects due to methanol crossover. At 
the anode side the presence of gas phase reduces the active cata- 
lyst areas, while at the cathode side the presence of liquid water 
reduces the active catalyst areas. The modeling results of polar- 
ization curves show good agreement with the experimental data. 


Compared with the single-phase model, the two-phase model 
shows better agreement with the experimental data, especially at 
high current densities. Further comparisons show that the single- 
phase model under-predicts the anode overpotential as well as 
the rate of methanol crossover, and the amounts of under predic- 
tion increase with current density. The modeling results of CO2 
mass fraction showed that the porosity of the anode diffusion 
layer played a very important role in the DMFC performance. 
With a low porosity, the produced carbon dioxide cannot be 
removed effectively from the catalyst layer, thus reducing the 
active catalyst surface area as well as blocking methanol from 
reaching the reaction zone. 
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